RELAXED EXCITED STATE GEOMETRIES 

IN A PEIERLS-HUBBARD MODEL FOR POLYENES 
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Using parameter values relevant to polyacetylene and its finite polyene analogs, we determine, 
via exact finite size diagonalization of a ID Peierls-Hubbard model, the relaxed geometry and optical 
absorption spectra of several states of interest in optical experiments (the singlet (V-Ag) and triplet 
(1 3 B U ) ground states, the singlet one {n l B u ) and two (2 1 A g ) photon gaps, and the triplet one (n 3 A g ) 
photon gap) in order to examine the question of how well this simple model can really capture the 
behavior of this simple material. 



INTRODUCTION 

Recent experimental and theoretical work have established for the whole class of exciting novel 
reduced dimensional materials, including the high-temperature superconducting copper oxides, "heavy- 
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fermion" systems, organic synthetic metals including conducting polymers, organic superconductors, 
and charge transfer salts, and halogen-bridged transition-metal chains, that it is essential to include 
both the electron-electron (e-e) and electron-phonon (e-p) interactions if we are to truly understand 
the properties of these materials. By carrying through detailed comparison of model calculations with 
real materials, we hope to achieve our modest goal: to gain an enhanced understanding of the nature 
of the competition between these two basic interactions of condensed matter physics. The hope of 
achieving our goal is greatest, of course, when the material and model are as simple as possible while 
still retaining the essence of the competition between the e-e and e-p interactions. Thus, studying a 
linear chain material of simple atomic constituents should be ideal. In fact, polyacetylene, or (CH)^, 
has previously been touted as the "hydrogen atom of the conducting polymers." 

No interpretation of dynamics of finite polyenes and conducting polymers can be considered com- 
plete unless it correctly captures the important role of the triplet and other excited states in these 
systems. Triplet states in (derivatives of) finite polyenes are known to be of considerable importance 
in the processes of vision [1] and of photosynthesis [2] , in that they control certain crucial relaxation 
paths after electronic excitation. In conducting polymers, triplets appear essential to the understand- 
ing of the nature of long-lived excited states and specifically of the spectra observed in photoinduced 
photoabsorption. From a theoretical perspective, knowledge of excited states and their properties 



can provide additional insight into the relative importance of electron-electron and electron phonon 
interactions and can help determine the parameters in the theoretical models of these materials. 

MODEL 

In the context of conducting polymers, the one-dimensional Peierls-extended Hubbard Hamilto- 
nian [3] provides a theoretical framework capable of treating both electron-phonon (e-p) and electron- 
electron (e-e) interactions of arbitrary strengths. We focus on the one-dimensional Peierls-extended 
Hubbard Hamiltonian (1-D PHH) in the context of irans-polyacetylene (trans- (CH) X ). The 1-D PHH 
is 
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where the general features of this Hamiltonian have been described elsewhere [3]. Our studies here 
are a continuation of the work described in Ref.s [4], [5], and [6] as to it's specific applicability to 
polyacetylene. We now can incorporate additional excited state information. 

Due to space limitations, we refer the reader to Ref.s [4] and [5] for our determination of polyene 
parameters using ground state data and to Ref. [6] for initial studies of excited geometries using 
conventional polyacetylene parameters. We found that, since the bandwidth is partly due to e-e and 
partly e-p interactions, the hopping integral is less than the conventional t$=2.h eV value, but the 
parameter ratios remain roughly the same. Our best fits were obtained with t$=\-\.h eV, a=3.5- 
4.5 eV/A, [7=3-6 eV, U=0-2 eV (with the fit roughly depending only on U-2V), and If =40-60 
eV/A 2 . A typical absorption using parameters in this range showing agreement with triplet and 
soliton absorptions was shown in Fig. 3 of Ref. [5] . 

EXCITED STATES 

Using the method for extracting excited states described in the appendix, we use the Lanczos 
method to determine the energies and wavefunctions of the ground state and several important excited 
states: the lowest triplet state (1 3 B U ), the second singlet "A s " state (2 1 A g ) which determines the 
two-photon gap, and the "optical states" which determine the singlet optical gap (n 1 !?^) and triplet 
optical gap (n 3 A g ) [7]. Initial studies using conventional polyacetylene parameters and a fixed uniformly 
dimerized geometry were reported in Ref. [6] . Initial relaxed studies of only the triplet geometry were 
also reported in that paper. We have now obtained relaxed geometry information for all these excited 
states. Both conventional polyacetylene parameters and our new set given above have roughly U/t =3 
and a<5/io=0.1, where S is the average dimerization in the ground state [8]. The relaxed geometries 
for this (dimensionless) parameter set are shown in Fig. 1 for a 12-site system using open boundary 
conditions. Rather than the raw distortion, we have plotted the distortion order parameter ±(—l) n 8 n 
(with the sign being chosen so that a short bond at the chain end is positive) so that the deviation from 
uniform distortion is more apparent. We also plot the total energy for each of these states at each of 
these geometries to show the expected luminescence shifts, as well as check that our numerical procedure 
converged to the lowest energy configuration and not a higher lying local minimum. Unfortunately, 
for this small an electron-phonon coupling, it is difficult to analyze the geometry dependence as it 
is so weak. Thus, following Ref. [9], we increase the electron-phonon coupling until a5/t =0.5. The 
results are shown in Fig. 2. Note that the singlet gap at the relaxed geometry of the optical state is 
approximately half its value at the ground state geometry. We also show in Fig. 3 the size dependence 
of the relaxed geometries for these same parameters. The state labeled "2 1 A g " is higher in energy than 
the state labeled "2 A*" on small systems (iV=8,10) but becomes lower in energy on larger systems 
(7V=12). Note that in the ground state geometry, the 2 1 A g lies roughly twice as far above the ground 



Fig. 1 (left). Relaxed geometry and total energy of the various ground and excited states discussed in 
the text for U/t =3 and a5/t =0.1. 

Fig. 2 (right). Same as Fig. 1, but for a5/t =0.5. For this electron-phonon coupling, a higher lying 
2 1 A g geometry can also be found. 

state as does the 1 3 B U state. This strongly hints that the 2 1 A g state is in some sense composed 
of two triplets, as suggested also by the work of Tavan and Schulten [10]. From Fig. 3 we see that 
the 1 3 B U state has the lattice distortion appropriate to two neutral solitons, and the 2 1 A g prefers a 4- 
soliton configuration at N=12. Does this mean that one should view the 2 1 A g state as composed of four 
solitons? Such an interpretation was proposed in the renormalization group studies of Hayden and Mele 
[11] but is explicitly contradicted in the configuration interaction studies [10] and (at least apparently) 
in combined experimental/theoretical studies on octatetraene [12]. Unfortunately, it is difficult to 
determine a definitive answer to this question (see Ref.s [3] and [6] for a more thorough discussion). At 
present, our data on the fully relaxed geometry and energy of the 2 1 A g state on the largest systems (N 
= 12,14) are not sufficient to determine the resolution of this question unambiguously. This difficulty 
arises in part because the 2 1 A g state is not the lowest singlet state of A g symmetry and thus must be 
extracted from the Lanczos data with considerable care. However, to date our Lanczos results are most 
consistent with the "4-soliton" viewpoint for long polyenes, and not inconsistent with the experiment 
showing a different geometry for octatetraene. 
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Fig. 3. Size dependence of the relaxed geometries for the parameters of Fig. 2. For iV=8,10 the 
configuration labeled 2 1 A* is preferred, but at N=12 the configuration labeled 2 1 A g becomes lower in 
energy. 

APPENDIX: NUMERICAL METHOD FOR OBTAINING EXCITED WAVEFUNC- 
TIONS 

Our "exact" diagonalizations were performed using a standard Lanczos algorithm which in essence 
involves expressing the Hamiltonian in a cleverly chosen basis, in which it is tridiagonal. Only the 
electronic part of the Hamiltonian was treated exactly; the adiabatic approximation was used for the 
phonons. Our formulation of the standard Lanczos procedure and method for extracting ground state 
relaxed geometry and optical absorption information, including our averaging technique for enhancing 
the finite-size results, has been previously described (Ref. [13]). 

While the standard Lanczos method is useful for projecting out ground states, we are typically 
interested in several states which are important for, e.g., optical experiments: the (ground state), 
2 1 A g (two-photon singlet gap), n x B u (singlet optical gap, \ X B U for small U), l 3 B u (triplet ground 
state), and n 3 A g (triplet gap, l 3 A g for small U). While the standard Lanczos method allows us to get 
the singlet and triplet (by working in the S z =l subspace) ground state wavefunctions and energies, 
it allows us to find only the energy of the one and two photon gap states, and not the wavefunctions 
of these excited states, and hence we cannot investigate, e.g., questions of geometry relaxation upon 
excitation and corresponding luminescence spectra. However, by slightly modifying the Hamiltonian 
and procedure, we can indeed obtain such information, though some care must be taken as occasionally 
the numerical procedure fails and so one must test that the resultant states are indeed eigenfunctions 
of the original Hamiltonian with the desired symmetry. Once they have been found, evaluation of 
correlation functions (such as the bond-charge, which determines the self-consistent relaxed geometry) 
proceeds as for the ground state. Though we describe here the method only for a single boundary 
condition, the generalization to our "boundary condition averaged" scheme, described in Ref. [13] is 
straightforward. 

The first step [14] is to note that the starting Hamiltonian, Eq. 1, commutes with several symmetry 
(parity) operators, especially if the distortion has some symmetry. Assuming there are equal numbers 



of up and down spin electrons and that the distortion is symmetric about some bond (or site), call it the 
n-th one, then we can add the following terms to the Hamiltonian without changing the eigenfunctions: 



X F F + X R R + X S S 2 



Where S is the total spin, R is the parity operator which inverts about the ra-th bond (site), and F 
flips every spin. When H = H + Hi is applied to a wavefunction, one finds 



Thus we can, e.g., preferentially select the singlet over the triplet ground state by setting X F to a large 
number and performing the standard Lanczos procedure. One could also apply a projection operator 
P which projects out a given total spin. However, this particular operator is (i) computationally 
expensive and (ii) tends to cause the numerical procedure to fail (through, e.g., loss of orthogonality). 
Thus we have not employed it. Further, raising or lowering the energy based on the total spin tends 
to favor zero (singlet) and maximum spin, rather than the singlet and triplet, and flipping the spin 
has the same singlet /triplet selectivity, thus we have only used the symmetries F and R to project out 
eigenstates, and tested the resultant state to insure the desired total spin. 

To project out the \ X B U is a simple matter of adjusting the A and finding the "ground" state 
(of H, first allowed 1 B U for H ) with the usual Lanczos procedure. However, for large U, the \ l B u 
is not the optical edge [7]. One can get around this problem by a slight modification of the Lanczos 
procedure. If we assume the ground state \4>q) is known, then after we have generated a given basis 
vector we also calculate and store (0o|J|^j). Once the storage limit of the computer is reached, 
instead of selecting the eigenvector with lowest eigenvalue, we select the one with lowest eigenvalue 
and non-zero (4>o\J\4>k) ■ (To do this we need only deal with the small matrix which converts a Lanczos 
basis number to a eigenvector number, thus this step does not add considerably to the time or size of 
the calculation). We iterate this procedure until the estimate of the n x B u energy stops changing, as for 
the ground state. Note that this (i) only increases the memory requirement by one state, and (ii) only 
adds one "expensive" step per Lanczos basis state. Of course, the statement "non-zero" above means 
in practice larger than some cutoff. If too small a cutoff is used, the l x B u , 2 1 A g , triplet ground state, 
or other lower lying state not of current interest, may be found instead of the desired optically allowed 
state. We have found a relatively large value of the cutoff (~ 10 -2 in units of the matrix element of the 
optical edge) is required, though this is not as large as it seems, as it is the square which is important 
for optical absorption. It also means that, e.g., for the triplet ground state geometry, it is not the band 
edge, but rather the soliton absorption, which is found (unless quite a large value for the cutoff is used, 
which can be troublesome, or an additional selection requirement that the energy be "near" the value 
for the uniform geometry, assumed to have been previously calculated, is added). 

At the end of this procedure for projecting out a given excited state, it is important to test 
that the resultant state (i) is an eigenstate of H , and (ii) has the desired symmetry and total spin, 
to insure against loss of numerical accuracy. (We have found using A values on the order of the 
energy separation of interest is usually successful.) Obviously, while straightforward, this excited state 
projection procedure cannot be used blindly, and it is important to test not only that the "answer" is 
an eigenstate with correct symmetry and spin (and non-zero (4>o\J\(j)k) if an optically allowed state is 
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E ( 1 A g ) + X F + X R 



Eq( 1 B u ) + Xp — Xr 
Eq( 3 B u ) — Xp — Xr + 2Xs 
E ( 3 A g )-X F + X R + 2X s 



desired) but that comprehensible trends as a function of system size or other parameter are followed, 
before one can have full confidence that one has indeed calculated what one set out to. 

When there are unequal numbers of up and down spin electrons or the distortion has no inver- 
sion symmetry, the symmetry operators F or R, respectively, cannot be used, though the procedure 
described above (with obvious modifications) can still be used to project out, e.g., the wavefunction 
of the first optically allowed state. Also it is clear that one can obtain further excited state informa- 
tion (e.g., second optically allowed state) using these and similar modifications of the basic Lanczos 
procedure, though at the cost of larger time and storage requirements and (potentially fatal) loss of 
numerical accuracy. We note that further excited state information, such as the many-body density 
of states, can also be obtained in a related but somewhat different fashion which essentially involves 
recursively generating the basis to fill in the necessary matrix element of the correlation function of 
interest and which does not significantly increase the storage requirements or decrease the numerical 
accuracy, though it obviously increases the time requirements [15]. 
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